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It is argued that continuum realisations of distributions of collisionless charged particles should 
accommodate a dynamically evolving number of electric currents even if the continuum is composed 
of only one species of particle, such as electrons. A model is proposed that self-consistently de- 
scribes the interaction of such a continuum and its electromagnetic field. It is formulated using a 
Lagrangian approach and employs a "folded" flow map to describe the bulk particle motion. An 
asymptotic perturbation scheme is developed to analyse ultra-relativistic multi-component current 
configurations. The model is fully relativistic and is formulated over Minkowski spacetime using 
intrinsic tensor field theory. 
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I. INTRODUCTION 

Progress in high-energy physics relies on accelerator 
designers contemplating charged particle beams of ever 
higher energies and intensities. As schemes for acceler- 
ating charged particles become more complex and am- 
bitious in their aims, it is apparent that some existing 
theoretical models are inadequate for a proper under- 
standing of new challenges. In many existing models, 
matter is represented using classical point particles and 
it is not clear how to unambiguously and consistently 
model their electromagnetic interaction. 

The nub of the problem is precisely how one should 
sensibly formulate the interaction of a classical point 
charge with its own electromagnetic field. The charge 
density of a point particle is singular and must therefore 
be handled carefully. Assumptions must be made about 
how the singular Coulombic stresses are compensated by 
non-electromagnetic stresses when calculating the par- 
ticle's self-force. This issue was addressed by Dirac pj 
nearly 70 years ago and led to the covariant Lorentz- 
Dirac equation for the trajectory of a point charge in an 
external electromagnetic field. However, unlike more fa- 
miliar equations of particle mechanics, the Lorentz-Dirac 
equation is a third- order ordinary differential equation 
in proper time for the particle's trajectory and has a 
number of unusual properties including self-acceleration 
and prc-accclcration. Unless special final conditions are 
adopted, it predicts that an isolated and free point charge 
with non-zero initial velocity will accelerate forever and 
if the particle is subjected to a sharp electromagnetic 
pulse it will begin to move before the pulse reaches it. 
Methods for evading such unpalatable conclusions in- 
volve iterating the Lorentz-Dirac equation in powers of 
its charge qo. Landau and Lifshitz |2| showed that trun- 
cating the resulting series to any finite order leads to a 
second-order evolution equation yielding particle trajec- 
tories with more reasonable properties. Although their 
argument yields a workable scheme, it is not at all clear 
that it may be generalised to a collection of accelerating 
high-energy charged particles in close proximity. Neglect- 



ing higher-order terms in qo may be suspect when the 
particle number density is sufficiently high. In conclu- 
sion, a number of ad-hoc choices must be made to obtain 
a sensible relativistic equation of motion for a collection 
of charged point particles starting from first principles. 
For a recent account of the issues concerning the deriva- 
tion of the Lorentz-Dirac equation see . 

Many of the above issues are due to the uneasy mar- 
riage of field and point particle concepts in classical elec- 
trodynamics. Classical point particles are a convenient 
idealisation and one might argue that quantum theory 
must be invoked to yield a palatable answer. However, 
the testy relationship between fields and point parti- 
cles resurfaces in quantum electrodynamics where diver- 
gences in "bare" quantities must be regulated prior to 
renormalisation to obtain physical results. 

Given the above complexities and reservations, an al- 
ternative approach has recently been developed [3| to 
analyse the ultra-relativistic dynamics of a collection of 
accelerating charged particles. The attitude adopted 
in Q and the present article is that models of matter 
based on classical relativistic field theory are more suit- 
able for high-energy beam dynamics than those employ- 
ing classical point particle notions. The model in 0] em- 
ploys a smooth relativistic field description of the elec- 
tromagnetic and matter content where the total energy 
and momentum of the electromagnetic and matter fields 
are conserved. Charged matter is modelled as a 4-vector 
field on spacetime whose trajectories describe the bulk 
particle motion. 

The partial differential equations governing the elec- 
tromagnetic and matter fields in this article are fully 
coupled and non-linear. Although exact solutions de- 
scribing highly symmetric configurations can be found, 
the system of equations is in general only tractable when 
subjected to an approximation scheme. Such a scheme, 
based on a covariant asymptotic expansion in a running 
parameter e > around the light-cone, was introduced 
in [3| and permits calculation of the detailed dynamics of 
ultra-relativistic charged particle beams in external elec- 
tromagnetic fields. In principle, one can calculate the 
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field quantities to any desired finite order in e. 

The present article focuses on some of the issues en- 
countered when analysing the model in Q. Specifically, 
the non-linearities in the field equations may lead to solu- 
tions whose charge density diverges despite being initially 
regular. Sections [H] and IIIII briefly review the charged 
continuum model in 3] and show that diverging solutions 
exist satisfying a substantial range of initial conditions. 
This behaviour is commonplace in many continuum mod- 
els in physics, such as in neutral gas dynamics and fluid 
dynamics, and is often ameliorated by including dissipa- 
tive processes. Section IIVI argues that dissipative pro- 
cesses will not prevent the formation of multiple currents 
in charged beams. SectionslVland lVll develoD and analyse 
a continuum model accommodating currents with a dy- 
namical number of components and section rVHI extends 
the ultra-relativistic analysis methods introduced in |3j 
to multi-component charged currents. 

II. SINGLE CURRENT CHARGED CONTINUA 

The model discussed in this section and in Q describes 
a collection of accelerating charged particles, with (rest) 
mass mo and charge go> as a dynamical continuum. The 
vector field V is the 4- velocity of the continuum on space- 
time and its integral curves describe the bulk motion of 
the collection of charged particles. The scalar field AT 
is the particle number density measured by a comoving 
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observer and p = j° 2 Af is called the reduced proper 

fcQ / / LQ 

charge density where e is the permittivity of the vac- 
uum and c is the speed of light in the vacuum. In what 
follows, units are chosen in which c = 1. 

The antisymmetric rank 2 covariant tensor field F (a 

2- form) encodes the electromagnetic field and the triple 
(V, p, F) satisfies the covariant Maxwell equations 3 

dF = 0, (1) 
d*F = -p-kV , (2) 

on Minkowski spacetime (A4i,g) where g is the met- 
ric tensor. In an inertial Cartesian coordinate system 
(t, x, y, z) in the laboratory frame 

g = —dt (g) dt + dx (8 dx + dy (g> dy + dz (g> dz, 

d is the exterior derivative, * is the Hodge map associ- 
ated with g and the 1-form V is defined by the property 
V{X) = g(V,X) for any vector field X. The field equa- 
tions for V are obtained using energy-momentum conser- 
vation drx = where the total stress-energy-momentum 

3- form tk is the sum of matter and electromagnetic con- 
tributions : 

T K = pg{V,K)*V + ^(i K FA*F-FAi K *F) 

and the vector field K is a, spacetime translation 011M4. 
Setting K to d t , d x , d y and d z in dru = yields the d t , 



d x , d y and d z components of the equation 

V V V= i v F, (3) 

where 

g(Y,V) = -l (4) 

with V the Levi-Civita connection on JH4 and iy the in- 
terior (contraction) operator on forms. The term iyF in 
@ is a continuum generalisation of the covariant expres- 
sion for the Lorentz force on a point charge (the charge 
to mass ratio qo/mo has been absorbed into the defini- 
tions of p and F) where the tangent to the point charge's 
proper-time parametrised worldline has been replaced by 
V. The charged matter drives the electromagnetic field 
through J5J and the electromagnetic field acts back on 
the matter through © conserving total energy and mo- 
mentum. 

Equations 111 II) are well-known in charged plasma 
physics and are often said to describe a "cold charged 
fluid". They have found application to accelerator 
physics in recent years and have proved useful for 
examining the stability of high intensity particle beams. 

III. DEVELOPMENT OF SINGULARITIES IN 
THE CHARGE DENSITY 

In Q highly symmetric exact solutions to (111 H describ- 
ing "walls of charge" were used to motivate a hierarchy 
of field equations for modelling ultra-relativistic charged 
particle beams. Exact solutions to (111 111 of the form 

F = £(t, z) dt A dz, 

v= . = =(d t + p(t,z)d z ) (5) 

y/l - p 2 {t,z) 

were sought, where pd z is the Newtonian velocity of the 
charge distribution as measured by the laboratory ob- 
server dt- Using l]l !■">!) it follows that 

dS = P W, (6) 
V V V = £#V, (7) 
g(V,V)=-l (8) 

where the projected metric g is 

g = —dt <g> dt + dz ® dz, 

# is the Hodge map associated with the volume 2-form 
#1 = dt A dz, V is the Levi-Civita connection of ~g 
and V(X) = g~(V, X) for any vector X on 2-dimensional 
Minkowski spacetime A^2 with metric g. Equation © 
implies 

V£ = (9) 

i.e. the electric field is constant along the integral curves 
of V and so, using (7J, the magnitude of the acceleration 
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Vy V is constant along the integral curves of V. There- 
fore, the continuum undergoes local hyperbolic motion 
and it is straightforward to solve to (jfilSJI in a comovin 
coordinate system (r, a) adapted to V, as shown in 
where V — d T and z = a on the initial hypersurface 
t = 0. For charge distributions initially at rest the Ja- 
cobian of the transformation between (r, a) and (t, z) is 
non-degenerate for all r and a [3j. Using it follows 
that, for all r > 0, p is well-behaved for charge distribu- 
tions at rest at r = 0. 

The purpose of this section is to demonstrate that more 
general initial conditions lead to divergences in the re- 
duced proper charge density p over finite time. In Q 
particular examples of p were generated using expressions 
for V and £ as functions of t = f(t, z) and a — a(t,z). 
For the present purposes it is more convenient to for- 
mulate an ordinary differential equation for p along V 
and examine properties of its solutions. The integrabil- 
ity condition 

d(p#V) = 
following from JSJ is written 

Vp = -p6 (10) 

where the scalar = ^~ 1 d#V is the divergence of V . 
Equations flT))) . JHJ and V9 = f(p,6,£) for some / 
forms a closed first-order ordinary differential system for 
(p, 9, £) along the integral curves of V. An explicit ex- 
pression for / is obtained below. 

Using the identity (see, for example, page 229 of Q) 

(C v g)(X,Y) =g(X,V Y V) +g(Y,V x V) (11) 

where X and Y are any vector fields on A^2 and Cy is 
the Lie derivative with respect to V, it follows that 



(C v V)(X) = (C v g)(V,X) 

= g(V,VxV)+g(X,V v V) 

and, using (JHJ), 

g(V,V x V) = ^7 x (g(V,V)) 
= 0, 



(12) 



(13) 



since V is metric-compatible. Therefore, (|12f) yields 

C V V = V V V. (14) 

Using Cartan's identity Q Cy = iyd + diy on forms it 
follows 



£ V V = i v dV 



(15) 



and so the following expressions for the relativistic accel- 
eration A = VyU of V are obtained: 



where (JTJl, (fHjl and ljTB|) have been used and A(X) = 
~g(A, X) for any vector X on M.2- Since there are no 
3-forms on 2-dimensional manifolds V A dV = and 



i v (V AdV) = 



(17) 



= (iyV)dV -V Ai v dV. 
Using jni, (Unj and i v V = g(V, V) = -1, it follows that 



dV = E#l. 



(18) 



Any frame field (Xq,Xi) and its dual coframe field 
(e^e 1 ) on M2 satisfy 



e a (X b ) = 51 

' ' \ if a = b 
if a j£ b 



5 a h 



(19) 



where the indices a, b run over 0, 1. The intrinsic curva- 
ture of V is zero so 

Vy V Xa V - V Xa Vy V - V[v,jr.] V = (20) 

where [V, X a ) is the Lie bracket of V and A a . The diver- 
gence 9 — # _1 d#U of V may be written J9J 



v- y = e a (y Xa v) 



hence 



V9 = Vy(V- U) 

= (Vye fl )(Vx a F)+e°(WVxJ) 

= (Vve")(Vx/) + VJ| e "(V [V jjy) 



(21) 



where (|20H and (|16|l have been used. The torsion of V 
vanishes 

V v X a -V Xa V-{V,X a }=0 
and so, using ifH?|) 

V v e a = -(e a (X7 v X b ))e b 



= -(e a (X7 Xb V)+e a ([V,X b }))e b . 
Using to eliminate Vye a in (|2T|) yields 
V6> = V ■ .4 - tr(WW) 



(22) 



(23) 



where W =_ e b {V Xa V)e a ® X b _and_tr(VFW) = 
e h (V Xa U)e a (Vx il U). The scalar tr(WW) is obtained 
using the following g-orthonormal frame field {Ao, Ai} 
adapted to V and its dual frame field {e°, e 1 }: 



An 



V, A 1 = #U, 



A = Vy V = i v dV = £#V 



(16) 



-V, 



(24) 
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where the vector field Xi = #V is defined by g(X, X\) = 
on any vector X onM 2 . Using (O, {TT} and 
(|2Ul it follows that 



tr(WW) = {g{X u W Xl V)f 



\({C v g){XuXi)) 2 



(25) 



and expressing g as 



3 = -e° <g> e° + e 1 ® e 1 

= -v <g> u + #u ® #u 



it follows 



and 



£ y g = - £ y V ® 1/ + F (g) £ y V 

+ £y#U <8> #U + #y <g) £y#U 

(£v^)(Jfi,Jf 1 ) = 2(£v#V)(X 1 ) 
= 26> 



(26) 



where JT£J, pi)l. = #- x d#V and Cartan's identity 
£y = iyd + diy on forms have been used. Equations 
(E3J), and J2U give 

U0 = V • .A - 6> 2 (27) 

and writing V • A as a differential form yields 



#- l d#i v dV 



(28) 



= -#- 1 rf(yA#rfV) 

= -#~\dv a #du) + #- x (y A d#dV). 

Thus, using ©, GEJ, and J2HJ it follows that is 

V6> = £ 2 + p-6l 2 . (29) 

Equations H29|) , (|1U|I and © are a closed system of dif- 
ferential equations for (p, 9, £) along V. 

Let r be any proper-time parametrised integral curve 
ofV: 

T:I->M 2 , 

A — > (t = T(A),« = Z(A)) 

where / is a subset of the real line R and 

r»d A = v. 

Hence, equations ©, JUJ and J2H) pulled back to R 
using T* are 



0, 



d£ T 
~dX 
dpr 
dX 

__£ r + pr -0 r 



-Prvr, 



(30) 

(31) 
(32) 



where the subscript T indicates pull-back using T* e.g. 
£ r (A) = (T*S)(X) = £(T(X),Z(X)). The general solution 
to JSJ is 

p r (A)=p cxp^-^ r (A')rfA'^ (33) 
where po = Pr(0) is a value of p on an initial hypersur- 
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rAf and the proper number density 



face. Since p 

TV > it follows pr > 0. Over intervals of A on which 
dpr/dX ^ (1371321 yields 



cfpr 

d6> r _ 

dpr -pr#i 



and so 



pr-0 2 ) 



2£ 2 



Pr^(pF 2 ^) = -^-2 
dpr Pr 



leading to the first integral 



9 2 = £'l + 2p r + K p 2 T 



(34) 



of 130I32|I where ko is a constant of integration deter- 
mined by the initial values #r(0) = #o> Pr(0) = Po and 
£r(0) - £b. 

According to l|34|l #p is a quadratic function in pr with 
So and Kq held fixed and the large A behaviour of pr 
and Or crucially depend on the sign of kq. Since £q is a 
positive constant and 9 r is positive (0r is real), it follows 
that if kq < then pr in l (3*3jl cannot be arbitrarily large. 



If no < then pr is bounded from above. 



(35) 



However, if ko > then no such bound on pr exists 
and, in principle, pr and 9r can attain arbitrarily large 
values. In fact, as will now be shown, if kq > then pr 
may diverge in finite proper time. 

Assume that the initial data satisfies kq > and 9q < 
0. Using (|34|1 to eliminate 6* 2 from the right-hand side of 
leads to 



do 



(36) 



and since pr > 0, kq > and 9q < it follows from l|36(l 
that 9y < 0. Therefore, using the negative root of l|34l) . 



l£% + 2p r + K pp, 

to eliminate 9r from the right-hand side of (|31|l 

dpr 
dX 



Pry£o+ 2pr + n p r 
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is obtained and pr asymptotes at proper time Aoo where 



Aoo = / 1 .> dp r . (37) 

J pr\/ £ o +2p r + KoPr 

pa 

The integrand is positive and (£q + 2pr) > so (pfTj) 
implies 

OO 

Aoo < / 7=r d Pr 

Po 

1 

i.e. Aqo is bounded from above. 

If K o > $o < at A = then pr diverges 
at proper time A = Aoo less than l/(poy^co). 

The constant kq may be obtained in terms of data on 
the spacelike hypersurface t = = T(0) where (i, z) is 
the laboratory coordinate system. Let p(0, z) be the ini- 
tial Newtonian velocity measured in the laboratory frame 
(see © for the definition of p) . Using QJ, 9 = V • V and 
it follows that 

e = d tl + d z { 1 p), (39) 
5 t7 + A*0*7 = (4°) 

where 7 = 1/ \J\ — p? . Using (|40|) to eliminate dtj from 
Pty yields 

= Ep + jd z p 

and so on the initial spacelike hypersurface t = 

6» = £oMo + 7o(9 2 Ai) (41) 

where p = (T»(0), (0,/*)„ = (T*(d z p))(0). Using gTJ) 
to eliminate 6*o = #r(0) in H34|) evaluated at A = gives 

K o = 4j (70(^)0 + 2£ po7o(^M)o - % - 2p ) ■ (42) 
Po \ 7o / 

Equation 142|l indicates that if (d z p)o = then k < 
and according to i|35[) pr does not diverge. This result 
agrees with the non-singular behaviour of the exact solu- 
tions to (|6I8J) satisfying p(0,z) = presented in y|; for 
Ko to be positive (d z p)o must be non-zero. 

The integral curves of V for the particular solution to 
equations (HilSIl with the initial conditions 

z oo 

2\J ^i-p(o, s y J v/i-Mo, s ) 2 / 

— OO 2 

p(0,z) = 0.01exp(-;z 2 ), 
p(0,z) = O.lsin(z) 



are shown in figure ^ Trajectories on which kq > and 
9q < are solid and trajectories on which kq < are 
dashed. Thus, according to l|38|l. p diverges at points 
along the solid trajectories, which meet at about t = 11. 
The details of the onset of the crossings are shown in 
ngure 121 with proper time A and z as axes. Evaluat- 
ing H37fl gives Aoo = 10.14 for the trajectory starting at 
(t, z) = (0,3). Any comoving coordinate system (r, a) 




z 



FIG. 1: The integral curves of V for a particular solution to 
equations 16181 . 

adapted to the trajectories, where a is constant on each 
trajectory, degenerates where the solid trajectories cross. 
The determinant of the Jacobian of the (r, a) — > (t, z) 
transformation vanishes at such points and p diverges. 



IV. DISCUSSION OF CROSSING 
TRAJECTORIES IN CHARGED PARTICLE 
BEAMS 

In fluid and gas dynamics, crossing trajectories are con- 
sidered to be a symptom of incomplete physics. All real 
fluids and gases are viscous to some extent and trajecto- 
ries may cross if their viscosity is neglected. For example, 
a compression wave in a hypothetical inviscid fluid may 
lead to a velocity field with crossing characteristics, i.e. 
the velocity becomes multi-valued. In reality this is not 
what happens; the velocity remains single- valued and sta- 
bilises to form a propagating shock. The fluid on either 
side of the shock is essentially inviscid but close to the 
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A 1 °- 




FIG. 2: The integral curves of V for a particular solution to 
equations 16181 . 



shock the second-order spatial derivatives of the velocity 
are so large that dissipation can no longer be neglected. 

However, it is far from clear that such arguments are 
relevant for a charged particle beam, which is physically 
very different from a normal fluid. Although (111 111 are 
often said to describe a "cold charged fluid" y| , this ter- 
minology is misleading. Microscopically, a normal fluid is 
a complicated system of neutral particles whose interac- 
tions are dominated by molecular collisions and possibly 
gravity but the dominant inter-particle forces in a beam 
of electrons are entirely electromagnetic in origin. At 
first sight, it seems that all of the necessary physics is 
contained in 111 II) . 

The velocity field V is a smoothed out representation of 
the particle motion. The fact that the trajectories cross 
does not mean that the particles are colliding; it merely 
indicates that a smooth field representation of the par- 
ticles has degenerated. Similarly, the 3-volume number 
density of a collection of particles may attain arbitrarily 
large values if the particles dynamically arrange them- 
selves into planar or linear configurations. 

Nevertheless, the electric field induced by V in ijll II) is 
inconsistent where trajectories of V cross. Equation @ 
indicates that the electric field in any wall-of-charge so- 
lution is constant along V. However, during the crossing 
the charge distribution "passes through itself" and the 
electric field has to change. The reason for this is easiest 
to appreciate by the following simple analogy. Consider 



the dynamics of a pair of positive sheet charges that are 
permitted to pass through each other. The sheets are 
arranged so that their normals and their electric fields lie 
along the z-axis. Each sheet has the same properties; its 
self-induced electric field is E at all points to its right and 
—E&t all points its left, where E > is constant, and its 
charge per unit area is Q = 2eoE. The sheets are labelled 
1 and 2 and sheet 1 lies to the left of sheet 2 initially. The 
electrostatic force acting on each sheet is constant and is 
due to the electric field of the other sheet. The force per 
unit area acting on sheet 1 is —QE and on sheet 2 is QE. 
Thus, the sheets will repel each other but their electric 
forces will not be enough to stop them meeting for suf- 
ficiently large opposing initial velocities. Let the initial 
velocities be so large that the sheets pass through each 
other. At all times before they meet, the force on sheet 1 
is — QE and the force on sheet 2 is QE. After they meet 
their roles have been exchanged; sheet 2 is acted on by 
— QE while sheet 1 is acted on by QE. Now consider a 
large number of sheets undergoing collective motions in 
which only some of the sheets pass through each other. 
A continuum realisation of this model is a dynamical set 
of component continua each with its own velocity field. 
The number of components evolves in time and depends 
on the history of the total continuum. 

The spacetime fields (V, p) satisfying 111 II) offer an Eu- 
lerian description of a single component continuum. Al- 
though a dynamical number of components can be simu- 
lated using more complicated Eulerian field theories [||, 
such approaches are restrictive because an upper bound 
must be placed on the number of anticipated compo- 
nents. In the following section a new Lagrangian model 
of a multi-component charged continuum is proposed in 
which the number of components is dynamical and free 
to attain any value. The essential idea is to describe the 
bulk particle motion using a flow map C from a body- 
time manifold into spacetime rather than inducing the 
motion from a velocity field V on spacetime. C may 
be described as "folding" a single electric current on the 
body-time manifold to give a multi-component current 
on spacetime. 



V. LAGRANGIAN DESCRIPTION OF 
MULTI-COMPONENT CHARGED CONTINUA 

Ingredients in the following Lagrangian description are 
an auxiliary 4-dimensional manifold B± called the body- 
time manifold and a map C from B4 into Minkowski 
spacetime M4. The body-time manifold B4 — K x B3 
where B3 is a 3-dimensional body manifold so each point 
P G B4 is also written P = (A,P) for A G K and PeB 3 . 
Each point P G B3 generates a curve Cp in spacetime 
M.4 where 

Cp{X) = C(\,P). 
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The map C is normalised so that A is the proper-time 
parameter of Cp_ for all P € S3: 



g(C,C) = -l 



(43) 



where C(P) = (C*d x )(P) is a vector at p = C(P) in M4. 

In general, C is a many-to-one map, i.e. there exists 
Pi and P 2 in £4 such that C{P X ) = C(P 2 ), and C is 
not required to be surjective. For any point p G A4± 
there may exist any number N(p) of real roots of the 
equation p = C(P). The map C describes "multi- valued 
velocities" because although C{P\) = C(P2), there is no 
reason why C(P\) should equal C(P2)- Thus, in general 
C cannot be identified with a vector field on ^4. The 
domain of C is B4 and C(P) is a vector at the point 
C(P); the map C is referred to as a vector field over C. 

The map C is defined to satisfy the Lorentz force equa- 
tion 



which vanishes at critical points. At each point p in M4 
the partial 4-current Jui (p) has the form 

J[i\ (p) = Q(Pm)C(P[{\ ) where Pfa = (p) (49) 
and p is the scalar field 

6 = TXT (5°) 



on 23. The system 1)431461) differs significantly from 0} 
because the number N(p) of elements of C -1 ({p}), 
and therefore the number of partial 4-currents in 14611 . 
depends on the spacetime point p. Numerically integrat- 
ing (|43I46|) involves computing N(p) at each time step. 



VI. MULTI-COMPONENT CHARGE 
CONFIGURATIONS 



V d C = i d F 



(44) 



where F is an electromagnetic field 2-form on M.4 and 
(V ( j,C')(P) and i ( j^F(p) are covectors at p = C(P). 

The maps V^C* and i^F are covector fields over C (i.e. 
1-forms over C). 

The set inverse C^ 1 of C at p includes the set of points 
in $4 for which p = C (P) and is written 



c- x ({ P }) 



if N(p) 



■,P[N( P )]} ifJV(p)>l 

: 



where the square brackets distinguish root labels from 
coordinate and frame labels. 

Each element of C~ 1 ({p}) gives rise to a partial electric 
4-current J[i}(p), where i = 1, 2, N(jp). The sum of 
partial currents is the total electric 4-current driving F 
in the Maxwell equations 



dF = 0, 



d*F 



N(p) ^_ 

* J [< 



(45) 
(46) 



i=0 



Regions of spacetime with different numbers of par- 
tial currents are distinguished by examining the pull-back 
C*(*l) of the spacetime volume 4- form *1 by C. Critical 
points in B4 are defined by the vanishing of the 4-form 
C*(*l). Their images under C are also said to be critical 
and lie in the interfaces between spacetime regions with 
different N(p). Specifying a non- vanishing closed 3-form 
J on #4 satisfying 



i dx J = 7 £ 9x J = 
leads to a scalar field A on 84 where 
AdAAj = C*(*l) 



(47) 



(48) 



Since (111 111 are equivalent to (|43I46|) when applied to a 
spacetime region with a single partial current, it follows 
that solutions to 111 II) and I|43I46II agree in such a region. 
This is illustrated in figures and 0] by collapsing "wall 
of charge" solutions to 1)11 II) and 143I46[) . The ansatze 
for F and V are 

F = £(t, z)dt A dz, 

V = coshx(i, z)d t + sinhx(i, z)d z 

and p depends only on (t, z) where (i, x, y, z) is an inertial 
Cartesian coordinate system with 

g = —dt ® dt + dx ® dx + dy (£> dy + dz ® dz. 

Similarly, C and g only depend on (A, a 1 ) where 
(a 1 , a 2 , a 3 ) is a coordinate system on B3. The initial 
conditions on £, \ and p are 



£(0,z) = ^( / l/M-nsh \ j(l)..sW/> •- / (/m>.sl)\ )(().. s) ( /.s- ). 



0.025 if -1.5 < 3 < 1.5 



J otherwise, 
%(0,z) = -1.2tanh(z) 



(51) 



with analogous initial conditions on C and g. 

In both cases a critical point develops at p cr it = 
(0,1.075). The integral curves of V in the single- 
component Eulerian model 1)114(1 exhibit crossings along 
a narrow corridor inside the forward light-cone of p cr it 
(see figure |3J). However, the solution to (|43I46(I shown in 
figure 01 contains a "fan" of three partial currents inside 
Pcrit's forward light-cone. The two models clearly yield 
dramatically different results. 

On the other hand, the two models agree outside the 
forward light-cone at p C rit- An argument for this is as 
follows: Integrate (111 H and I|43I46|I using a time slic- 
ing adapted to a field of synchronised inertial observers 



moving along the z-axis in the laboratory frame with 
positive constant velocity arbitrarily close to the speed 
of light. As before, the solutions agree up to a constant 
proper time surface containing the point p C rit- Part of 
this proper time hypersurface almost coincides with the 
z > subset of Merit's forward light cone. Now integrate 
the equations using a time-slicing adapted to a field of 
synchronised inertial observers moving along the z-axis 
in the negative direction at almost the speed of light. The 
result agrees almost up to the z < subset of p C rit's for- 
ward light-cone. It follows that the solutions to ijll II) and 
(|43I46[) agree at points outside the forward light-cone at 

Pcrit- 

For further comparison the reduced proper charge den- 
sities for 111 II) and (|43I46() are shown in figures [5] and 




z 



FIG. 3: The trajectories of the particular solution to Q 
with the initial conditions I51D . Trajectories cross in a nar- 
row region inside the forward light-cone of the critical point 
p cr it = (0,1.075). 



VII. ULTRA- RELATIVISTIC 
APPROXIMATION SCHEME 

Ultra-relativistic solutions to (111 II) may be obtained by 
promoting (V, p, F) to a 1-parameter family (V £ , p £ , F £ ) 
in e where the e dependences are motivated by exact 
"wall of charge" solutions The equations obtained 



S 




z 



FIG. 4: The trajectories of the particular solution to 11431 
1461 1 with the initial conditions 15111 . The forward light-cone of 
the critical point p cr it = (0, 1.075) is also shown. The region 
outside the "fan" emanating from p cr it has only one partial 
current whereas the region inside the "fan" contains three 
partial currents. 



by equating orders in e lead to a self-consistent hierar- 
chical method for approximating solutions to (111 If) . The 
virtue of the scheme is that the e dependences conspire to 
produce an infinite tower of equations that are partially 
coupled and are generally easier to solve than the fully 
coupled system (111 H . 

A similar approach based on l|43l4(j[) will now be out- 
lined. Let C' £ be a 1-parameter family of maps from B4 
into M4 such that 

g(C £ ,C £ ) = -l 

where C £ = C £ d\. Let X — es and introduce the map C £ 
where 

C e {s,P) = C £ (es,F) (52) 
at all non-critical points in the domain of C £ . Thus 

C £ = -^ds = Jc £ ' (53) 

where C £> = C £ d s and so 

g(C £ \C £ ') = ~e 2 . (54) 



3 1.5 



3 1.5 



FIG. 5: The reduced proper charge density of the particular 
solution to (111 It with the initial conditions 1511 . 



and the Lorentz force equation for C e becomes 

V C e,C*' = ei C z<F e . (55) 
In Q a dependence for V E on e of the form 



E ^ 

n=-l 



(56) 



was motivated by "wall of charge" solutions to (11111) . 
For multi-current configurations on Minkowski spacetime 
M4, the flow map C is the dependent variable and it is 
natural to exploit the afline structure of M.4 and postu- 
late an analogous series for C e . Let (x a ) be an inertial 
Cartesian coordinate system adapted to the laboratory 
frame on M.4 where a, b = 0, 1, 2, 3 and 



where 



9 = Vabdx a ® dx b 



-1 ifa = 6 = 
Vab = <i 1 if a = 6 7^ 
if a ^ b. 



In the rest of this article the map C s is also regarded as 
the 4-component column vector 




FIG. 6: The partial reduced proper charge densities of the 
particular solution to 1431461 with the initial conditions (15 II . 
The region outside the "fan" has only one partial current 
whereas the region inside the "fan" contains three partial cur- 
rents. 



where (C ea ) are the (x a ) components of C £ 
The e expansions of C sa (Q), where Q 
chosen as 



00 

C-(Q) = ^>"OQ) 




(s,P), are 



(57) 



Motivated by the corresponding expression in the Eule- 
rian single-current formulation [3j , the 1-parameter fam- 
ily F 6 of electromagnetic 2-forms is chosen as 



F e {p) 



E 

n=-l 



e n F n (p) 



(58) 



at any point p in M4. Note that (|58fl is independent of 
the e expansion of C e in l|57|) ; F £ is a 1-parameter family 
of 2-forms on Minkowski spacetime M4 and l|58|) is not 
the e expansion of the 2-form F e (C £ (Q)) on B4. 

One way to minimise the complexity of the ensuing 
calculation is to adapt a coordinate system on B3 to J. 
Let (C 1 , C 2 , <^ 3 ) be a coordinate system on B3 and let J7i23 
be the corresponding component of J: 

J = Jm^ 1 Ad£ 2 Ad£ 3 . 



(59) 



Eliminating g in Ju 
and l|59fl yields 



in favour of J and C using (|48I5U|) 
i7m(-P[i] 



det(DC)(P w ; 



C(P W ) 
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where P w = CT^{p) = (A, £\ £ 2 , £ 3 ) and DC is the Jaco- 

bian of C a (X, £ x , £ 2 , £ 3 ). By definition, J123 is indepen- 
dent of A (see so (£ 1 ,£ 2 ,£ 3 ) may be chosen so that 
|Ji 23 | = 1. Hence 



and 



J[i\(p) 



1 



C{P {i] ) 



|det(DC)(P H )| 
1 

\{det{T>C)oC^){p)\ 



and so 

•TO 



|(det(DC-)oC f 7 1 )(p)| 



r'oC^ 1 )' 1 ^) (60) 



where (l")2t and (f 531) have been used, DC e is the Jacobian 
of e EO (s,f\£ 2 ,£ 3 ) and DC e = ±DC e . 

Let F h eQ = r] ac F§ c where F e = \F e ah dx a A efe fc and 
ff c f]cb = &t where 



81 



1 if a = b 
if a ^ 6. 



The inertial coordinate representations of the Lorentz 
force equation H55f) for C £ and the normalisation con- 
straint (|54ll are 



eai nebi 



C ea " = e(F £a o C E )C E( ", ryahC ea, C 
where F £ is a solution to the Maxwell equations 

N(p) ^ 

dF e = 0, d * F e = - * J [{\ 

1=0 



(61) 



(62) 



and Jg, is given in JfJUJl. 

Inserting (|S7jl and into (|B*T]) and (|B^|l and equating 
equal order terms in e induces a hierarchy of equations 
for successive approximations to (C £ ,F £ ). The first six 
steps in the hierarchy are: 

• Adopt an external electromagnetic field i.e. a 
solution to the source-free Maxwell equations 

dF_i=0, d*F_i=0. 

• Solve 

Cg"(Q) = (F- 1 ) b a (p)C b '(Q) 

for Cg subject to 

r lab Q / (Q)C b o / (Q) = 

where p = Cq(Q) and (F_i)b a is data obtained in 
the previous step. 

• The 2-form Fo is a solution to the Maxwell equa- 
tions 

N(p) 

dF = 0, d * F = - ^2 * J Wo 

i=0 



J$ ]0 (p) = \det(VCZl)(p)\C$ (Q) (63) 



where DCj^p) is the Jacobian of C,^ at p = (x a ) 
and Q = C^ip). 

• The first order correction C\ to the map Co is ob- 
tained from the linear equation 

Cf = (F^) b a C h ^ + (f- 1 ) 6 a )C CfCg' + (F ) b a Co" 

subject to 

= 

where maps on 2?4 are implicitly evaluated at Q and 
maps on M.4 are evaluated at p = Cq(Q). Indices 
following a comma indicate partial differentiation 
with respect to the corresponding coordinates so 
(F-i)b a ,c = ■^{F^ 1 ) b a . 

• The 2-form F\ is a solution to the Maxwell equa- 
tions 

N(p) 

dFi =0, d-kF% = — 

i=0 

and 

J^=|dct(DC W0 )| -^(DC^DCx)^ 
+ {tr(DC W0 D 2 C ) T DC H0 C 1 }C 

- a 

+ c;-Dc Dc H0 d 

where maps on M.4 are implicitly evaluated at p 
and maps on #4 are implicitly evaluated at Q — -- 
CZ Q (p)- Inside the square brackets DCi is the ma- 
trix of derivatives of Cf and DC is the matrix of 
derivatives of C' a and both should be regarded as 
linear maps from the (£° = A,^ 1 ,^ 2 ,^ 3 ) compo- 
nents of vectors on B4 to the (x a ) components of 
vectors on M.4. Similarly, inside the square brack- 
ets DCj~j should be regarded as a linear mapping 
from the (x a ) components of vectors on M.4 to the 
(£ a ) components of vectors on B4. The column vec- 
tor tr(DC [ 7 ]0 D 2 C ) is 



/ tr(DC H0 D9 A C ) \ 
tr(DC^ D9 e Co) 
tr(DC^ D%Co) 

V tr ( DC Ho D % C o) J 



where Dc\Co is the matrix of partial derivatives of 
cVCo whereC = A,£\£ 2 ,£ 3 . 
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• The second order correction C 2 to C is a solution to 
the equation 

Cf'(Q) = (F^) b a C b 2 ' + {F^) b \ c CiC b ^ 

1 ^~ I T? \ a pcpdpbl 1 / p \ a pcpbl 
+ -^\r-l)b xd<--i<--i<-o +l-^-lJb ,ct-2<-"0 

+ (f )6 o Cf + (^o) b Q xCi c C 6 ' + (F 1 ) 6 °C 6 ' 
subject to 

27 7afc C a 'C 2 b ' + ? 7afc C 1 a 'C 1 b ' = -l 
where p = Co(Q). 

Note that C and Cq + eC^ are lightlike and C„ + eC( +£ 2 C 2 
is the leading order timelike approximation to C e . 

In single-current regions the reduced proper charge 
density is ^/— g( J^, J^) on spacetime. Using (|o7|) . ljoTj|) 
and the normalisation condition (|61fl on C £/ it follows 
that J— g(J^, J^) converges to as e tends to 0. The 

reduced proper charge density diverges on interfaces be- 
tween single-current and multi-current regions and the 
above approximation scheme is valid arbitrarily close to 
such interfaces. 

VIII. CONCLUSION 



normal fluids are viscous to some extent and ameliorate 
the problem by forming a shock. On the other hand, a 
beam of electrons is very different from a normal fluid 
and it was argued that there is no physical reason why 
the trajectories modelling an electron beam cannot cross. 
A Lagrangian theory permitting crossing trajectories was 
presented and its "wall of charge" solutions were exam- 
ined and compared with solutions to the original Eulerian 
field system. Finally, an approximation scheme was de- 
veloped to analyse ultra-relativistic charge configurations 
of the Lagrangian system. 

The Lagrangian theory discussed here features an TV- 
component electric current where N is dynamically de- 
termined and has a point-wise dependence on spacetime. 
It is possible that configurations with N = 1 initially 
may evolve into highly complicated "turbulent" configu- 
rations where N is arbitrarily large. Further work in this 
context may be found in Q. 

The theory developed here has immediate application 
to high-energy accelerator physics where ultra-relativistic 
motion is ubiquitous. The general approach is also valid 
in systems where electromagnetic interactions dominate 
over collisional processes, such as in laser-driven plasma 
wakefield accelerators where the transition between 
single and multiple component electron currents ("wave 
breaking" ) Q has recently received much attention. 



Continuum models of charged particle beams include 
the interaction of matter with its own electromagnetic 
field and avoid peculiar phenomena evident in point- 
particle descriptions of self-interacting charge, such as 
self- acceleration and pre-acceleration. A field-theoretic 
realisation of a collection of classical electrons is the 
"cold" charged continuum. It was shown that the veloc- 
ity field of the continuum may possess crossing trajecto- 
ries and in this case the Eulerian theory is inconsistent. 
Such behaviour is not uncommon in physics; compres- 
sion waves in inviscid fluids also develop crossing trajec- 
tories. However, inviscid fluids are an idealisation; all 
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